subroutine denomt1(t1, r1, ie, sym)
use mod_iop
use mod_size
use mod_orbit
use mod_ioff 
implicit none
integer, intent(in) :: sym
real*8, intent(inout) :: t1(*), r1(*), ie(*)
!
if(prec==8.and.igpu==0) then 
   call denomt1_d(t1,r1,ie,sym)
elseif(prec==4.and.igpu==0) then 
   call denomt1_s(t1,r1,ie,sym)
elseif(prec==8.and.igpu==1) then 
   call denomt1_dg(t1,r1,ie,sym)
elseif(prec==4.and.igpu==1) then 
   call denomt1_sg(t1,r1,ie,sym)
endif
!
return 
end
!
subroutine denomt2(t2, r2, ie, sym)
use mod_iop
use mod_size
use mod_orbit
use mod_ioff 
implicit none
integer, intent(in) :: sym
real*8, intent(inout) :: t2(*), r2(*), ie(*)
if(prec==8.and.igpu==0) then 
   call denomt2_d(t2,r2,ie,sym)
elseif(prec==4.and.igpu==0) then 
   call denomt2_s(t2,r2,ie,sym)
elseif(prec==8.and.igpu==1) then 
   call denomt2_dg(t2,r2,ie,sym)
elseif(prec==4.and.igpu==1) then 
   call denomt2_sg(t2,r2,ie,sym)
endif
!
return
end
